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We have performed a single crystal neutron scattering experiment on Ceo.7Lao.3Be to micro- 
scopically investigate the order parameter of phase IV. Below the phase transition temperature 
1.5 K of phase IV, weak but distinct superlattice reflections at the scattering vector k = (4, 5, 5) 
(h, I — odd number) have been observed for the first time by neutron scattering. The intensity 
of the superlattice reflections is stronger for high scattering vectors, which is quite different 
from the usual magnetic form factor of magnetic dipoles. This result directly evidences that 
the order parameter of phase IV has a complex magnetization density, consistent with the 
recent experimental and theoretical prediction in which the order parameter is the magnetic 
octupoles T 13 with Ts symmetry of the point group Oh ■ Neutron scattering experiments using 
short wavelength neutrons, as done in this study, could become a general method to study the 
high-rank multipoles in f electron systems. 
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The importance of high-rank multipolar degrees of 
freedom of f electrons in strongly correlated electron 
systems has recently been widely recognized. A typi- 
cal dense Kondo compound CeBg with a simple cubic 
crystal structure of space group Pm3m is a well-known 
example where the importance was clarified experimen- 
tally and theoretically. It shows the following two suc- 
cessive phase transitions: The first is from the paramag- 
netic phase (phase I) to the antiferro-quadrupolar order- 
ing phase characterized by the wave vector &q = [5,3,5] 
at 3.3 K (phase II), followed by the second transition 
to the antiferromagnetic ordering phase with a complex 
magnetic structure characterized by the four nonequiva- 
lent wave vectors at 2.3 K (phase III), 1 where the fifteen 
multipoles in the F§ quartet crystal-field ground state of 
CeB6 play an important role for these orderings. 2,3 By 
doping La into the Ce site in this system, a new phase 
called phase IV appears below Try = 1.7 K and 1.5 K 
in Ce^Lai-^Be for x = 0.75 and 0.70, respectively. 4-11 
From the initial discovery of phase IV, 4-6 its character- 
istic magnetic phase diagram and anomalous bulk prop- 
erties, which show the isotropic cusp of magnetization 
at Try and the strong elastic softening of C44 within 
phase IV, suggested that the ordering of phase IV must 
be different from any quadrupolar ordering. Since the T$ 
quartet has three types of octupoles - T xyz , T a , and T 13 
— in addition to dipoles and quadrupoles, it was argued 
that the magnetic octupoles are a possible candidate for 
the order parameter. Following these studies, the magne- 
tization under uniaxial pressure, thermal expansion, and 
elastic constant measurements supported that the order 
parameter of phase IV is T' 3 octupoles with T 5 symmetry 
of the point group G\; 8 ~ n it is consistent with the exis- 
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fence of internal mag netic fields detected by NMR 12 and 
/iSR 13 ' 14 as well as a theoretical model. 15 Furthermore, 
evidence of the antiferro-octupolar ordering of T° with 
the same wave vector as A;q has recently been reported 
by the resonant X-ray scattering experiment and its de- 
tailed analysis, 16, 17 although further studies are neces- 
sary for understanding the overall nature of phase IV in 
CexLai-zBg. 8 

In principle, such time-reversal-symmetry-breaking 
high-rank multipoles can be detected by neutron scatter- 
ing, because a neutron interacts with electrons through 
magnetic interactions. 19 ' 20 Moreover, this probe has the 
advantage of being able to get direct information about 
the magnetization density of high-rank multipoles from 
the magnetic form factor. Notwithstanding this expecta- 
tion, no significant magnetic Bragg peaks were observed 
for phase IV within experimental accuracy in our pre- 
vious neutron scattering experiments below M = ^j^- 
— 0.4A -1 , where k, 9, and A are the scattering vector, 
Bragg angle, and neutron wavelength, respectively. 2 1-23 
Although the contribution from magnetic octupoles to 
the neutron scattering cross section may be small, it must 
be finite and the scattering intensity in the high k vec- 
tor region is expected to be stronger than that in the low 
k vector, 20 ' 24 because magnetic octupoles must have a 
complex magnetization density with no spatially uniform 
magnetization. Furthermore, a recent theoretical calcu- 
lation of the magnetic form factor of octupoles predicts 
a detectable scattering intensity magnitude. 24 

In this letter, we report the recent results of the neu- 
tron scattering experiment on a Ceo.7Lao.3B6 single crys- 
tal, focusing attention on superlattice reflections in the 
high scattering vector region. We have succeeded in de- 
tecting weak but distinct superlattice reflections from 
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phase IV by neutron scattering for the first time. The 
k dependence of the magnetic form factor in the super- 
lattice spots directly evidences that the order parameter 
of phase IV has a complex structure of magnetization 
density, consistent with the theoretical and experimen- 
tal prediction that the order parameter is the magnetic 
octupoles T° with symmetry of Oh- 

A large single crystal of Ceo.7Lao.3Bg was grown by 
the floating zone method, using 99.52% enriched 11 B to 
avoid the large neutron absorption due to 10 B. The bulk 
properties of the single crystal were checked by electrical 
resistivity and magnetization measurements. The sam- 
ple is cylindrical in shape with a diameter and length 
of 4.4 mm and 14 mm, respectively. The cylinder axis is 
nearly parallel to the [010] direction. The neutron scat- 
tering experiment was performed on the thermal neu- 
tron triple-axis spectrometer TOPAN (6G) at the JRR- 
3M reactor in the Japan Atomic Energy Agency. The 
sample was mounted in the mixing chamber of a 3 He- 
4 He dilution refrigerator with a superconducting magnet. 
Magnetic fields were applied along the [1,1,0] direction, 
normal to the (h,h,l) scattering plane. Incident neutrons 
with the short wavelength A = 1.4133A were selected 
by a pyrolytic graphite (PG) monochromator in order to 
search for superlattice reflections in the high scattering 
vector region. The triple-axis mode was used with the 
collimation open-60'-60'-60' and double PG filters to get 
a better signal-to-noise ratio. In this experimental setup, 
the mosaicity of the sample is 0.36° full width at half 
maximum (FWHM), reflecting the good quality of the 
single crystal. 

For determining of the magnetic form factor, some cor- 
rections are needed. The Lorentz factor and absorption 
factor corrections were made for the observed nuclear 
and superlattice reflections. The former was represented 
by 1/ sin 29 in the present geometry of scans. The latter 
was numerically calculated by the Fortran program by 
taking the approximate shape of the sample into account. 
The change in the intensity due to the absorption correc- 
tion is less than 16%. To obtain the absolute value of the 
magnetic form factor, information about the normaliza- 
tion factor between the integrated intensity of fundamen- 
tal nuclear Bragg reflections and the intensity calculated 
from the nuclear structure factors is also needed. For cal- 
culating the nuclear structure factors, the most reliable 
site parameter of B determined by the previous powder 
neutron experiment on Ceo.75Lao.25Bg 23 was used. How- 
ever, the integrated intensity of several nuclear Bragg 
reflections is not proportional to the calculated value, es- 
pecially for strong Bragg reflections. This deviation may 
be caused by the unavoidable extinction effects, 25 which 
strongly influence the Bragg intensity in experiments us- 
ing a large single crystal with a small mosaicity as used 
in this study. It is difficult to correct for the influence of 
the extinction effects on the nuclear Bragg intensity in 
this case. Thus, the normalization factor was estimated 
by using the intensity of three relatively weak Bragg re- 
flections (2,2,0), (1,1,0), and (0,0,1). 

At the lowest temperature 0.25 K under a zero mag- 
netic field in phase IV, we have observed weak superlat- 
tice reflections characterized by the wave vector 
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Fig. 1. (Color online) Scattering patterns of rocking curves at the 
scattering vector re = in phase I (at temperature T = 

2 K and magnetic field B = T), phase III (T = 0.25 K, B = 
1 T), and phase IV (T = 0.25 K, B = T) in Ceo.7Lao.3Be. The 
magnetic field is applied along the [1,1,0] direction. 
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Fig. 2. (Color online) Difference diffraction patterns between 
0.25 K and 2 K under a zero magnetic field at re = 
along the [1,1,1] direction in Ceo.7Lao.3Bg. The lines are Gaus- 
sian fits. 



for the first time. This wave vector is the same as de- 
termined by the recent resonant X-ray scattering exper- 
iment. 16 Figure 1 shows an example of the scattering 
patterns of rocking curves in phases I, III, and IV. The 
enhancement of the scattering intensity in phase IV has 
been clearly observed. The peaks seen in phases I and 
III are the extrinsic higher-order contamination from the 
strong (5,5,5) nuclear Bragg reflection. It should be noted 
that the observed distinct enhancement of intensity does 
not come from the extinction effects, which produce an 
increase of about 1% in the nuclear Bragg reflections in 
phase IV, 21 because the extinction effects do not affect 
such a weak signal and the enhancement of the inten- 
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Fig. 3. (Color online) (a) Temperature dependence of the inten- 
sity at the peak positions of re = (5,5,5), (§>§>§), an d (f if if) 
in a zero magnetic field, (b) magnetic field dependence along 
the [1,1,0] direction of the intensity at the peak position of k = 
(f ,f if ) at 0.25 K, where the background intensity is subtracted. 
The arrows indicate the phase boundary of phase IV reported by 
the bulk measurements. 



sity of the superlattice reflections is at least an order of 
magnitude larger than that of the nuclear Bragg reflec- 
tions. Therefore, we believe that the observed peaks are 
an intrinsic signal from the order parameter of phase IV. 
Figure 2 shows the difference diffraction pattern between 
0.25 K and 2 K under a zero magnetic field at K = 
( |,|,2) along the [1,1,1] direction with threefold symme- 
try. The lines are Gaussian fits. In the fitting procedure, 
the width of the Gaussian profile function is fixed to that 
of the fundamental Bragg reflections 0.36° (FWHM). We 
also observed several superlattice reflections along other 
directions of k, such as (f,f but reliable data with 
sufficient statistical accuracy have not been obtained be- 
cause of the very weak signal. Taking the Lorentz factor 
correction into account, we note that the integrated in- 
tensity at k = (§,§,§) and (§,§,§) is stronger than that 
at k — (5,0) 5); this is different from the usual magnetic 
form factor, as discussed later. 

To confirm that the observed signal comes entirely 
from the order parameter of phase IV, we measured 
the temperature and magnetic field dependences of the 
peak intensity of the superlattice reflections, as shown 
in Figs. 3(a) and 3(b). In the temperature dependence 
under a zero magnetic field, the peaks at k = (5,5,5), 
(§,§,§), and (§,§,§) develop below T w = 1.4 K with de- 
creasing temperatures. In the magnetic field dependence 
along the [1,1,0] direction at 0.25 K, the peak intensity 
at k = (5,5,5) seems to disappear above about 1 T, 
which agrees with the phase boundary of phase IV under 
magnetic fields along the [1,1,0] direction. 7 These results 
clearly show that the observed weak superlattice peaks 
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Fig. 4. (Color online) Magnetic form factor at the superlattice 
spots along the [1,1,1] (open circles) and [1,1,1] (filled circles) 
directions in phase IV of Ceo.7Lao.3B6. 



come from the order parameter of phase IV. 

In general, the neutron can be scattered by time- 
reversal-symmetry-breaking multipoles due to an elec- 
tron with spin s and momentum p through the mag- 
netic interaction, because the neutron has the magnetic 
dipole moment. The cross section of the single electron 
is proportional to |(A|Q^|A}| 2 , where |A) is a state of the 
electron. The scattering operator Q± is defined by 



Qi^ = exptin ■ r) < k x (s x k) : — -k x p 

' h\n\ 



(1) 



where r denotes the position of the electron and k — 
k/\k\. 20 Q± is related to the Fourier transform of the 
magnetization density M(n) of the electron as follows: 



2Q±Mb = — « x (M(k) x k). 



(2) 



The expectation value of 2Q± is the magnetic form fac- 
tor. Since the magnetization density of octupoles is com- 
pletely different from that of dipoles, as seen in eq. (2), 
their magnetic form factors show a qualitatively differ- 
ent behavior in the k dependence. In the dipolar case, 
the form factor decreases with increasing n, such as that 
reported in the pure system CeBg. 25 On the other hand, 
in the octupolar case, the form factor is zero at k = 
and has a maximum at a finite k, reflecting its complex 
magnetization density with no spatially uniform magne- 
tization. The characteristic behavior of the form factor 
of octupoles is given by the detailed theoretical calcu- 
lation of 2|(A|Qj_|A)|, taking the orbital contribution of 
the second term of eq. (1) correctly, in which the atomic 
wave function diagonal for the magnetic octuples T' 3 in 
the Tg, quartet is assumed. 24 Therefore, we can distin- 
guish whether the observed superlattice reflections orig- 
inate from dipoles or octupoles by the k dependence of 
the magnetic form factor. 

Figure 4 shows the magnetic form factor at the super- 
lattice spots along the [1,1,1] and [1,1,1] directions ob- 
tained by using the integrated intensity of the observed 
superlattice reflections and by making the corrections. 
The absolute value of the form factor is slightly lesser 
than 0.1/iB- This magnitude especially at the smallest 
scattering vector k = (5,5,5) with = O.lA -1 , which 
is expected to be the largest scattering intensity in the 
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usual magnetic neutron scattering, is quite small; it is 
comparable to the magnitude of tiny magnetic dipole 
moments reported in some heavy electron systems. This 
smallness must be the reason why the magnetic reflec- 
tions were not previously detected even by the highest 
intensity powder neutron diffraction. 23 The form factor 
in Fig. 4 is clearly strong for high scattering vectors. 
This cannot be explained by the usual antiferromagnetic 
ordering even by considering any conceivable magnetic 
structure as well as any domain distribution because the 
data in Fig. 4 are the form factors at the superlattice 
spots along the same direction. Therefore, this unusual 
k dependence of the form factor directly evidences that 
the order parameter has a magnetization density differ- 
ent from ordinary dipole ordcrings. This result qualita- 
tively agrees with the theoretical calculation considering 
an average of four domains of the order parameter T 13 . 
Furthermore, from the selection rule of the cross sec- 
tion based on the symmetry classification of octupolar 
scattering for the three possible octupoles T xyz , T a , and 
T^, 24 the present result can rule out the possibility of 
T xyz with T 2 symmetry because the superlattice reflec- 
tions along the [1,1,1] direction with threefold symmetry 
have been observed. T a with T4 symmetry is also unlikely 
to explain the observed k dependence because magnetic 
dipoles with the same T4 symmetry as T a are expected 
to be mixed. Therefore, the present result strongly indi- 
cates that the order parameter of phase IV is the mag- 
netic octuples T 13 with T 5 symmetry. For a quantitative 
comparison between the experimental form factor and 
the theoretical calculation including the anisotropy of the 
form factor, detailed corrections of the extinction effects 
and data at more superlattice reflection points as well as 
more statistical accuracy are necessary. 

It still remains possible that such a k dependence of in- 
tensity might arise from lattice distortions, since the scat- 
tering intensity due to lattice distortions also increases 
with increasing k in proportion to the square of k. This 
possibility cannot be completely ruled out only on the 
basis of the present neutron scattering experiment, but 
it should be noted that no superlattice reflection due to 
lattice distortions has been observed by X-ray scattering, 
which is the most powerful probe for detecting lattice dis- 
tortions. This fact as well as the existence of the internal 
fields detected by NMR 12 and ^SR 13 ' 14 support that the 
observed superlattice reflections are magnetic. To defini- 
tively confirm this, we will perform a polarized neutron 
scattering experiment soon. 

In summary, weak but distinct superlattice reflections 
at k = (h, I = odd number) from phase IV of 

Ceo.7Lao.3B6 have been observed for the first time by 
elastic neutron scattering using a large single crystal. 
The intensity of the superlattice reflections is stronger 
for high scattering vectors. This unusual n dependence 
of the intensity evidences that the order parameter of 
phase IV has a complex magnetization density, consis- 
tent with the recent experimental and theoretical pre- 
diction in which the order parameter is the magnetic oc- 
tupoles T 13 with T 5 symmetry of Oh ■ Neutron scattering 
experiments using short wavelength neutrons, as done 
in this study, could become a general method to study 



the high-rank multipolcs in f electron systems. Further 
neutron scattering measurements, including experiments 
under uniaxial stress, are planned to clarify the detailed 
nature of the order parameter of phase IV. 
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